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Abstract 

We report results of studies of inclusive B — > nK + ir~ and B — > rfK + ir~ decays. Charged conju- 
gates are implied throughout this paper. These are obtained from a data sample containing 386 mil- 
lion BB pairs, collected at the T(45) resonance, with the Belle detector at the KEKB asymmetric 
energy e + e~ collider. The branching fraction of inclusive B° — ► rjK + iT^ and B° — > 7]TT + ir^ are mea- 
sured to be B(B° -► nK + ir-)= (31.7 ± 1.9±|;§) x 1(T 6 and B(B° -► rpr + ir-)=(6.2±\-$}%$) x KT 6 , 
where the first error is statistical and the second systematic. The decays B° — > a^X + , where 
X + = K + , 7r + were searched for and no significant signals found. Upper limits of B(B° — ► 
a K + )< 1.6 x 1(T 6 and B(B° -> a "vr + )< 2.8 x 10 -6 at 90% C.L. are obtained. Here the notation 
B(B° — > Oq X + ) indicates the product of branching fractions for B° — ► a l + and a$ — ► 777T - . 

PACS numbers: 
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INTRODUCTION 

Recently, observations of large branching fractions of B mesons to three-body charmless 
hadronic systems have been reported by the B factory experiments In the mesonic 

decays B + — ► K + -k + -k~ and B + — > K + K + K~ , the broad K + K~ mass spectrum above 1.5 
GeV/c 2 in B + — ► K + K~K + suggests a large non-resonant _B + — ► K + K~K + contribution. 
In the baryonic decay B + — > ppK + , the mass spectrum cannot be explained by a simple 
phase-space model. A baryonic form factor model[5] or an additional unknown resonance 
around 2 GeV/c 2 are both possible explanations. These studies of three-body decays have 
already provided new information on the mechanism of B meson decay. Further, they 
suggest opportunities to search for previously unknown decays. Here we report results on B 
meson decays to t]K + it~ and 7]7t + tt~ based on a data sample that contains 386 million BB 
pairs, collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 GeV on 
8 GeV) collider Q. KEKB operates at the T(4S) resonance (y/s = 10.58 GeV) with a peak 
luminosity that has exceeded 1.5 x 10 34 cm~ 2 s -1 . 

There are several interesting quasi- two-body decays such as B° — > a^K + and B° — > n + 
included in this study. The observation of these modes will provide information both on B 
meson decays to scalar mesons and on the nature of the %. 



APPARATUS AND DATA SET 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere p}. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
sample of 152 million BB pairs (Set 1), while a 1.5 cm beampipe, a 4-layer silicon detector 
and a small-cell inner drift chamber were used to record the remaining 234 million BB pairs 
(Set 11)0. 

For Monte Carlo (MC) studies, samples of signal, generic b— >c decays and charmless 
rare B decays are generated with the EVTGEN|9J] event generator. Continuum MC events 
from the process e + e~ — > 7* — > qq are generated with the JETSET[l(j generator. The 
GEANT3 [ll| package is used for detector simulation. Signal features are studied with one 
hundred thousand MC events of each decay mode. Background B decays are studied with 
706 million generic BB events and "rare B" MC corresponding to 21 times luminosity, where 
the latter are charmless hadronic and radiative B decays with branching fractions taken from 
the Particle Data Group 2 • 



EVENT SELECTION AND RECONSTRUCTION 

Hadronic events are selected based on the charged track multiplicity and the total visible 
energy sum, which gives an efficiency greater than 99% for generic BB events. All primary 
charged tracks are required to satisfy track quality cuts based on their impact parameters 
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relative to the run- dependent interaction point (IP). The deviation from the IP is required 
to be within ±0.1 cm in the transverse direction and ± 2 cm in the longitudinal direction, 
where the z (longitudinal) axis is taken as opposite to the positron beam direction. Due 
to the detector response, a polar angle cut ( —0.866 < cos 9 < 0.956 ) is applied to all 
charged tracks. Only tracks not identified as muons or electrons are used in this analysis. 
Particle identification (PhD) is based on C K / (C n + L K ) information, where £k(tt) stands 
for the likelihood for charged kaons (pions). PID cuts are applied to all the charged par- 
ticles in this analysis. Unless otherwise stated, these are £k/(£tt + £k) > 0.6 for kaons 
and £k/(£tt + -^k) < 0.4 for pions. The PID efficiencies are 85% for kaons and 89% for 
pions, while the fake rates are 8% for pions faking kaons and 11% for kaons faking pions, 
respectively. Photon energies are required to be greater than 50 MeV (100 MeV) within the 
acceptance of the barrel (endcap) ECL. 

Candidate rj mesons are reconstructed through 1] — > 77 and i] —> tt + 7t~7t°. For B re- 
construction, the momenta of rj candidates are recalculated by applying the r\ mass con- 
straint in a vertex-mass constrained fit. For rj — > 77 decays, candidate rfs are selected with 
| cos#*| < 0.90, where 9* is the angle between the photon direction in the rj rest frame and 
the f] momentum in the lab frame, to suppress the soft photon combinatorial background 
and B — > K*^ feed-across. The rj mass regions are 0.500 GeV/c 2 - 0.575 GeV/c 2 for rj — > 77, 
and 0.535 GeV/c 2 - 0.560 GeV/c 2 for r] -> tt+tt-tt . 

K*° mesons are reconstructed from K*° — > K + tt~. Candidate ag mesons are recon- 
structed from Oq — > f]' K ^ and are required to have masses within 150 MeV/c 2 of the nominal 
value. For — > t]tt~ , the modulus of the helicity(| cos^^-^) must be less than 0.8, where 

^hei(a _ ) i s the angle between the ir + and B° in the ag rest frame. 

Note that, in the inclusive -B — > r]K + n~ channel, we have applied a D° — > K + n~ (1.841 
GeV/c 2 < Mktt < 1.888 GeV/c 2 ) veto and in the inclusive B° — > ?77r + 7r~ channel, we have 
applied D° -> 7r + 7r-(1.80 GeV/c 2 < M wn < 1.94 GeV/c 2 ) and D~ -> ^"(1.80 GeV/c 2 < 
M m < 1.94 GeV/c 2 ) vetos. 

B meson candidates are identified using the beam-energy constrained mass Mt, c = 

V^bcam - I-PbI 2 and the energy difference AE = E B - E heSLm , where -E be am = 5.29 GeV, 
and (Pb, Eb) is the four-momentum of the B candidate in the T(4S) rest frame. We impose 
a 2-D box cut on M bc and AE with M bc > 5.2 GeV/c 2 and \AE\ < 0.3 GeV for further 
analysis. A signal region with M bc > 5.27 GeV/c 2 and -0.1 GeV< AE < 0.08 GeV is se- 
lected to plot the projections of fits. A sideband region is defined with M bc < 5.26 GeV/c 2 
inside the box region. 

At most one candidate per event in each mode is required. The best candidate is chosen 
based on the sum of x 2 of the r] vertex-mass constrained fit and the x 2 of the Kir(in{) vertex 
fit. 

The dominant background for the three body B decay events comes from e + e~ — > qq 
continuum events, where q = u,d,s or c. In order to reduce this background, several shape 
variables are chosen to distinguish spherical BB events from jet-like continuum events. Five 
modified Fox- Wolfram moments [l3j and a measure of the momentum transverse to the 
event thrust axis (S±) 14] are combined into a Fisher discriminant 0]. The Probability 
Density Functions (PDFs) for this discriminant and cos8b, where 9b is the angle between 
the B flight direction and the beam direction in the T(4S') rest frame, are obtained using 
events in the signal and sideband regions from MC simulations of signal and qq events. The 
displacement along the beam direction between the signal B vertex and that of the other B, 
Az, also provides separation. For B events, the average value of Az is approximately 200 fim, 
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while continuum events have a common vertex. Additional discrimination is provided by the 
6-flavor tagging algorithm [16[ developed for time-dependent analysis at Belle. The flavor 
tagging procedure yields two outputs: q (= ± 1), which indicates the flavor of the tagging B, 
and r, which ranges from to 1, and is a measure of the confidence of the tag. Events with 
high values of r are well-tagged and are less likely to originate from continuum production. 
Thus, the quantity q ■ r can be used to discriminate against continuum events. The PDFs 
derived from the Fisher discriminant, the cos 9b distributions and the Az distributions are 
multiplied to form a likelihood ratio TZ— C s / (C s + £ q q), where C s (L q g) is the product of the 
signal (qq) probability densities. We achieve continuum background suppression by imposing 
q ■ r-dependent TZ requirements, based on a study of the signal significance (Ns / a/ Ns + Nb) 
using a MC sample, where Ns and Nb are signal and background yields, respectively. The 
effect of the TZ cut is studied by comparing the cut efficiency on reconstructed B + — > D°% + 
events in data and MC, for different values of TZ. A systematic error of ~2% is obtained for 
the TZ cut. 

The resolution of the signal M^ c width (o"A/ bc ) is verified using data and MC samples of 
reconstructed B + — > r]'K + , r]' — > 777r + 7r~ events. Our MC underestimates <JM hc by 8.98^5;5q % 
and 7.561q 49 % f° r the se ^ I and set II data samples, respectively. An inclusive r( sample, 
where rj' — > 7]7t + tt~ , is used to check the AE resolution. The ratio of data to MC AE widths 
in the 77 -> 77 and 77 -> n+n-jr modes are 1.03 ± 0.02 (1.05 ± 0.05) and 1.07 ± 0.03 (1.14 
± 0.05) for the set I(set II) data sample. 



ANALYSIS PROCEDURE 

Signal yields are obtained using an extended unbinned maximum likelihood (2-D ML) fit 
to the Mb c and AE distributions in the M^ c — AE box region after the TZ cut is applied. 
For N input candidates, the likelihood is defined as 

f -(N S +N C + N bh +Nra) N 

L(N S , N C ) = T7| IIIWs. + N oP<k + N bl P bh + N ra P rai ], (1) 

where P$ v Pd, Pbli an d Pra, are the two-dimensional probability densities for event i to 
be the signal, continuum, charm containing (b —>■ c) B decay backgrounds and charmless B 
decay backgrounds in the variables M^ c and AE, respectively. Poisson fluctuations for Ns 
and Nq are considered in this type of likelihood. For backgrounds other than continuum, 
N bb - and N ra are obtained from MC samples. Uncertainties in the MC PDFs are included in 
the systematics study. The continuum, b — > c and charmless B decay background PDFs are 
all obtained from the respective MC samples. 

The Mbc and AE shapes from continuum MC events are modeled by an ARGUS func- 
tion jl7] with a fixed end point at 5.29 GeV/c 2 and by a 2nd order Chebyshev polynomial, 
respectively. The shapes of signal, BB, and other rare charmless B decays are modeled by 
2D smooth functions. The AE distribution is found to be asymmetric, with a tail on the 
lower side due to 7 interactions with material in the front of the calorimeter and shower 
leakage out of the back side of the crystals. As a result, the AE resolution and the tail dis- 
tribution strongly depend on the 77 energy. In the inclusive B° r]K + TT~ and B° — > i]7i + ii^ 
studies, the 77 energy distribution for the signal events is not known apriori so we divide the 
data into three samples: P v < 1 GeV/c, 1 GeV/c < P v < 2 GeV/c and P v > 2 GeV/c. 

For decays with more than one sub-decay process, the final average results are obtained 
by fitting the sub-decay modes simultaneously with the expected efficiencies included in the 
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fit and with the branching fraction as the common output. This is equivalent to minimizing 
the sum of x 2 — —2 hi(£) as a function of the branching fraction over all considered sub- 



decay channels. The statistical significance (E) of the signal is defined as J — 21n(Lo/L max ), 
where L and L max denote the likelihood values for zero signal events, and the fit number 
of signal events, respectively. The 90% C.L. upper limit is calculated by finding x max such 
that 

jo L ( x ) dx _ qn% (?) 



MEASUREMENT OF BRANCHING FRACTIONS 



We correct our signal MC efficiency for several observed differences between it and data. 
Differences in the PID efficiencies are corrected using a D* + — ► D°tt + , D° — > K~tt + control 
sample. A study using D* partial reconstruction is used to correct the tracking efficiency 
To correct the 7r° reconstruction efficiency, a high momentum inclusive r\ sample is used 
where the ratios of reconstruced rj — > 7r 7r 7r° to reconstruced rj — > 77 in data and MC are 
compared. The MC simulation for low energy photons is further tested by comparing the r\ 
helicity distribution for data with MC predictions. The ratio of the single ir° reconstruction 
efficiencies in data as compared to MC is 0.924 with a conservative systematic error of 
3%. The TZ cut efficiency correction is determined using B + — > D°n + decays. The high- 
momentum r) sample study is also used to correct the efficiency of 77 reconstruction and mass 
cuts. The correction factor for the mass cut of 77 — > 77(7^71-°) is 0.990 ± 0.001 (0.993 ± 0.003) 
All examined efficiencies show fairly good agreement between data and MC samples. The 
tracking, PID, tt° and rj reconstruction efficiency systematic uncertainties are also obtained 
from the above studies. 



Inclusive B° — > r]K + ir and B° — > r]7r + ir 

Signal yields obtained from the two dimensional M^ c and AE extended unbinned max- 
imum likelihood fits are shown in Table UJ The backgrounds from generic BB decays and 
other rare B backgrounds are considered in the fit. For the mode B° —>■ rjK + n~ , the largest 
component of the rare decay background is B — > r]'K. For the sample of B° — ► r]7r + 7r _ , most 
of the rare B events come from B~ — ► rjn^, B~ — > p^r] and B° —>■ rjK*°. We use different 
2D smooth functions as signal PDFs in the different regions. Fig. ^ and Fig. El show the 
AE and M^ c projections for the entire rj momentum range and the three sub-ranges. The 
fit results are summarized in Table d 

The yield of B° — > K*r] accounts for about half of the inclusive B° — > t]K + 7T~ yield. The 
other half is not well understood so the default signal efficiency, e sig , is based on a model 
with 75% B° K*rj and 25% phase space. The difference between e sig for this model and a 
model with 50% B° — > K*r) and 50% phase space is included in the systematic error. In the 
case of B° — > r]7T + TT~, there is no obvious enhancement in the two body mass spectra and 
the yields in different rj momentum ranges are similar to the expectation from phase space 
B° — > r]7T + 7T~ . Therefore we use a phase space model of B° — > 7]7t + tt~ to determine e S i g for 
inclusive B° — > 7]7t + tt~ . 
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FIG. 1: Projections onto M\, c and AE from the extended unbinned 2D ML fit for inclusive B — > 
r]K + 7r~. The fit results with different r/ momenta ranges are shown in the different plots. (a),(b) 
: P v < 1 GeV/c, (c),(d) : 1 GeV/c < P v < 2 GeV/c, (e),(f): P v > 2 GeV/c, (g),(h) : whole P v 
region. Gray dashed lines show the continuum qq contribution; Turquoise dotted lines show the 
generic BB background, yellow shaded parts show the other rare B events; the red curves show 
the signal component; the blue curves show the sum of all above contributions; and the histograms 
show the data distribution. 



TABLE I: The fit yields from the two dimensional M^ c and AE extended unbinned maximum 
likelihood fit for inclusive B° — > t]K + tt~ and B° — > 7/7r + 7r~. 



r\ momentum(GeV/c) 



B°- 




, v - 


■* 77 


B°- 




, v - 


-> 3tt 


B°- 


-> T]TT + Tr~ , 


7] - 


♦ 77 


B°- 




rj - 


* 3vr 



P„ < 1 1 < P v < 2 



P n > 2 whole 



25.9 t\H 122.9 tlg° 7 409.7 ± 28.7 558.6 t 

3.8 ± 7.7 26.6 ig'jj 111.4 j^jj 141.8 1 

7.0™ 44.3™ 32.0™ 83.3 ± 2 

11.1 ±f;f -3.0 ±1;7 4.4 1^ 12.4 ±\ 



Two body mass spectra of r]K + n and r]7r + ir in the final state 

The Dalitz distributions are densely populated close to the edges, as can be seen in Fig. 01 
indicating possible resonances in our sample. We investigate these by examining the 1-D 
projections of the Dalitz plot. We divide each 2-body mass distribution into 100 MeV/c 2 
bins. A 2-D fit is then applied in each bin to determine the B yield. To prevent cross talk 
from resonances in the Kir region, we require the Kir mass to be larger than 2 GeV/c 2 
when showing the Kr\ and nr/ mass distributions |l8j . For the -B — > r]K + 7i~ final state, the 
signal PDF is based on a mixture of B° — > t]K* and B° — > i]K + 7T~ MC samples. Fig. 0] 
shows the signal yields obtained from the 2-D fits as functions of the 2-body masses. There 
is an obvious enhancement from the i^*°(892) resonance and an excess in the Kir mass 
region between 1.4 GeV/c 2 and 1.7 GeV/c 2 which may be due to 1^(1430) or #£(1430). 




FIG. 2: Projections onto M\, c and AE from the extended unbinned 2D ML fit for inclusive B — > 
rjir + TT^ . The fit results with different rj momenta ranges are shown in the different plots. (a),(b) 
: P v < 1 GeV/c, (c),(d): 1 GeV/c < P v < 2 GeV/c, (e),(f) : P v > 2 GeV/c, (g),(h) : whole P v 
region. Components are the same as in Figure 1. 
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FIG. 3: Dalitz plots of (a)B° — ► rjK + ir and (b)B — ► r/7r + 7r candidates from the i? signal region. 
We have applied the D veto in both plots. 



For B° — > i]7i + 7i^, the signal PDF is based on phase space MC. Fig. 0] shows the results. 
In the 7r + 7r~ mass distribution, there is a small enhancement at low mass, which may be 
due to either the p° or / (980). In the 7]ti ± mass distribution, there is a small excess in the 
1 GeV/c 2 region that may be from the Oq. This is discussed in greater detail in the next 
section. 
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FIG. 4: The left hand plots are signal yields from the 2-D fit as a function of (a)M K + w - , (c)M K + v 
and (e)M 7T ~ r] with the B° — > r]K + ir~ final state. The right hand plots are signal yields from the 
2-D fit as a function of (b)M 7r + 7r -, (d)-/^/^- and (f)M r?7r + with the B° — ► r]ir + ir~ final state. 
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In the B° — > Oq ag — > ^tt~ mode, the background from other rare B decay modes 
is included only in the rj — ► 77 subdecay mode, otherwise it is ignored due to the negligible 
contribution predicted from MC. For the rj — ► 7r + 7r~7r° subdecay mode, the generic BB 
background is neglected for the same reason. Figure Efa) and fc>) show the M^ c and AE 
projections from the extended unbinned 2D ML fit for B° — > t]K + tc~, with M vn - in the 
an mass region. From the fit, the yields for the rj — > 77 and the 77 — > 7r + 7r~7r° subdecay 



^3.8 
-2.55 



respectively. These yields are the sum of B 



-R+ and 



modes are 17.3l|;o an d "4.7 
non- resonant -B — > r]K + 7r~. 

In the B° —>■ ag tt + , — > 77tt~, 77 — ► 77 mode, we use smoothed 2-D histograms to model 
the BB and rare .B backgrounds. The fractions of these backgrounds are fixed from the MC. 
In the 77 — ► 7r + 7r _ 7r° mode, we consider only the b — > c background and neglect backgrounds 
from other rare B decays. Figure 0(c) and Efd) show the M^ c and AS projections from 
extended unbinned the 2D ML fit for B° — > i]7T + 7T~ , with M^- in the ag mass region. From 
the fit, the yields for the 77 — > 77 and the 77 — > 7r + 7r~7r° subdecay modes are 22.6 Ig^ 1 and -0.4 
:|g, respectively. Again note that the yields are the sum of B° — > %7r + 
B°-^r]7i + - 
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SYSTEMATIC ERROR 

Systematic errors may arise from the efficiency corrections and from the fitting process. 
The main sources of uncertainties in the efficiency corrections are from the reconstruction of 
low-momentum charged tracks, low-energy photon finding, and TZ cut efficiency, each at the 
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level of a few percent. The systematic errors include contributions of 2% for the 1Z cut, 1% 
per reconstructed charged particle for tracking, 0.5% per charged particle for PID, 3% for 
7T° reconstruction, and 3% for r\ reconstruction. The fitting systematic errors are estimated 
by varying the fitting function PDF variables by ± 1 a from the measured values. The 
variations in the fitted yields are then quadratically summed to get the total fit systematic 
uncertainty. The contributions to the systematic error are given in Table ITT1 

The intrinsic width of is 57 MeV/c 2 for the MC samples of B° — > a^K + and B° — > 
a^7r + used to calculate the efficiencies. If instead we assume the intrinsic width of clq to be 
100 MeV/c 2 [l2|], the signal efficiency decreases by 9.7%. This uncertainty is considered in 
the upper limit calculation. 

DISCUSSION AND CONCLUSION 

In summary, we have searched for inclusive charmless hadronic B° — > 7]K + tt~ and B° — > 
?77r + 7r~ decays and have observed the branching fractions listed below, we also give upper 
limits at 90% C.L. for the inclusive B° — ► ?77r + 7r~ : 

B(B -> r]K + ir-) = (31.7 ± 1.9+H) x 10~ 6 , 
B[B -> 777T+7T-) = (6.2±i;|±g;|) x 10- 6 (< 11.9 x lO" 6 ). 

We do not find significant signals for B° — > clqK + or B° — > tt + . The 90% C.L. upper 
limits on the respective branching fractions are: 

B(B° -> a K + ) < 1.6 x 10~ 6 , 
B(B° -> a 7T+) < 2.8 x 10~ 6 . 

The notation — > a X + ) indicates the product of branching fractions for B° — > a^X + 
and ap — > i]7T~ , where X + is a A" + or rr + . All the results are listed in Table ITTT1 

TABLE II: Summary of systematic errors for inclusive B — > rjK + 7r~ and inclusive B — > rjir^ir" . 
The systematics from 77 and 7r° reconstruction are cr^ and cr^o respectively; cr^ is the systematic 
from PID; &n b - is the systematic from the number of BB pairs in our data; a s i g is the systematic 
due to uncertainty in the signal PDF; a fa contains systematics from other fit component PDFs; oj r 
is the systematic from charged tracking; and or is from the 1Z cut. All are expressed as percentages 
(%)• 



decay mode 


<7tr <7R <?r, 






& sig 


c/it 


Sum 


B° -» ^jK+tt- 


2.0 2.1 3.0 


0.0 


1.5 


1.0 


4.0 


+3.3 
-4.6 


+6.9 
-7.6 




4.1 2.1 3.0 


3.0 


2.5 


1.0 


4.2 


+4.8 

-8.2 


+9.4 
-11.5 


B° -> r]K+Tr- 


2.5 2.1 3.0 


0.3 


1.6 


1.0 


4.1 


+3.4 
-5.4 


+7.1 

-8.2 


B° — > 77 77 7T + 7r~ 


2.0 2.0 3.0 


0.0 


1.2 


1.0 


0.0 


+5.2 
-4.8 


+6.8 
-6.5 




4.2 2.0 3.0 


3.0 


2.5 


1.0 


0.0 


+18.6 
-13.2 


+19.8 
-14.8 


B° -» ?77r + 7r" 


2.6 2.0 3.0 


0.9 


1.4 


1.0 


0.0 


+11.6 

-8.3 


+12.5 
-9.6 
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TABLE III: Fitting significance(£) and branching fractions(S) from the extended unbinned AE- 
Mb c 2-D ML fits. The significance of Og X + includes possible non-resonant events in the selected 
mass region. 



Mode 




£ 


£(io- 6 ) 


inclusive B° — 


-> f]K + TT 


24.0 


31.7 ± 1.9 tH 


B° -> a K + 




2.0 


(< 1-6) 


inclusive B° — 




3.6 


6.2 ±tg±85«11.9) 


B° -> a^TT + 




1.4 


(<2.8) 




FIG. 5: Projections onto Mb c and AE from the extended unbinned 2D ML fits for B — > Oq K + 
are shown in (a) and (b). The projectoions for B° — > 7r+ are shown in (c) and (d). Components 
are the same as in Figure 1. 
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